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Abstract 



We calculate the isospin symmetry violating effects to the reactions pp dvr"'" 
and np dir^ arising from the different hadron masses and from the Coulomb 
interaction between the positive pion and the deuteron. These effects are large 
I enough in the cross section and analyzing power Ay that they should be taken into 

' account in comparisons of accurate experiments in different charge channels. 

<N : 
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O I Isospin symmetry has been used with great success to relate the experimentally avail- 

^ I able pion production reactions 

pp — > dn^ (1) 

and 

np d-K^. (2) 

Another charge channel is nn dn^ , but this is impracticable. Further, positive (or 
negative) pion absorption on the deuteron 

Tc^d pp (3) 

can be related to (1) by detailed balance. Because the isospin one component of the 
np initial state has the probability |, then according to this symmetry the cross section 
for the reaction (2) is simply one half of the pp initiated reaction (1). Relative angular 
distributions and spin observables of both reactions should be the same. 

However, the isospin symmetry is not exact. Since the masses of the initial state 
nucleons as well as the masses of the charged and neutral pions are different, the = 
reaction (2) has a 5.9 MeV lower threshold (CMS) than the = 1 reaction (1). This 
has a potential to be a significant effect as compared with errors of modern accurate 
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experiments. Some of this difference can be accounted for by different phase space factors 
in the cross sections. In relative observables these factors cancel. 

After the removal of the trivial phase space factors, nevertheless, there remain some 
charge dependent effects in the transition amplitudes themselves. These have three basic 
origins: 

1) Kinematic difference: The relation of the initial momentum to the final momentum 
is not exactly the same for different reactions. Even for the same final momentum (the 
more relevant in threshold reactions with a large negative Q- value), the initial momentum 
in the baryon wave function is slightly different. The latter enters the transition matrix 
via the NN wave function in a way dependent on the partial wave and can, in principle, 
influence also the relative observables. 

2) The Coulomb interaction enters the channel (1) in both the initial and the final 
state, while it is absent in (2). Presumably for the high energy nucleons it can be 
neglected, but it can cause a significant correction in the civr"'" final state in the threshold 
region. One may safely neglect the magnetic interaction (although it may be relatively 
important in the small charge symmetry breaking in np c/vr"). 

3) Charge dependent strong force: A more interesting possibility is charge dependence 
in dynamics, i.e. the interaction involved. The lighter and longer-ranged 7r° exchange 
enters differently in the pp and np initial states. This gives rise to an isotensor force both 
in NN scattering and in the A^A^ —>■ AN transition potential 0] as well as in pion 
s-wave rescattering. Another source for a charge dependent dynamic force can be meson 
mixings tt^ ^ r] and p° ^ uj. However, the latter violate also charge symmetry and are 
presumably smaller. 

The aim of this Letter is to study the first two contributions to charge dependence in 
pion production. The motivation is an urgent need for an estimate because of experimen- 
tal reasons. The accuracy of experiments on np —>■ drc^ is aiming at a level where even 
the possible variation from the above sources may be significant in comparing different 
charge channels 0. Precision measurements of np scattering at intermediate energies 
often use pion production and isospin symmetry in normalizing the luminosity and the 
above differences could be a source of a systematic error. This in turn may have impli- 
cations e.g. to the determination of the (charged) pion coupling to nucleons, a topic of 
much controversy recently 0] . Also there is progress in threshold experiments, where the 
Coulomb corrections are most relevant 0, ^, 0, |^. In other reactions, such as pp —>■ ppir^ 
its inclusion has been essential 0. Therefore it is necessary, if not interesting, to get a 
quantitative estimate of charge dependence in the A^A^ — *• d7r reactions. Although these 
violations of the isospin symmetry may be small, to our knowledge there have not been 
any systematic calculations of these to really raise comparisons of experiments above ad 
hoc reliance on charge independence. 

We shall discuss the charge dependent effects in two parts, showing the change in the 
observables most likely to be compared (the differential cross section and the analyzing 
power Ay). First we calculate the kinematic effect due to the change of the initial 
state making all comparisons for the same final pion momentum in the CMS. Then 
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we proceed to incorporate the final state Coulomb corrections in different approximate 
ways. The first one is simply to apply the Coulomb penetration factors |C/p, and in the 
second one we replace the plane wave pion in the matrix elements by Coulomb scattering 
wave functions. Finally, the effect of the finite charge distribution is also estimated. In 
principle, one should consider also the polarizability of the deuteron. However, we end up 
with the conclusion that this is most likely a small effect, much smaller than the already 
small ones that we include. Technically including the polarizability would be also much 
harder. 

A detailed discussion of the model used can be found in Ref. In addition to 

the direct production from the initial A^A^ states, we include by the coupled-channels 
method p-wave pion rescattering through the A(1232) isobar (dominant above 350 MeV 
in p-wave production from ^D2) and also s-wave rescattering (dominant at threshold). 
The Reid soft core wave function is used for the deuteron. This model gives rather 
good predictions to the cross sections and spin observables through the A region, but 



overestimates threshold production by about 50 % |TT . 

In the naive quark model the mass difference between the A^A^ channels and the AA^ 
channels does not depend on the charge 0, because the same six quarks are involved in 
both. This means that the kinetic energy required for the excitation of a AA^ intermediate 
state is the same for the np and pp initial states. Therefore no difference in the transition 
A^A^ — > AA^ arises from the baryon masses at the same incident nucleon kinetic energy. It 
was also pointed out in Ref. that in lowest order the internal Coulomb effects should 
have a similar effect to the total masses in all channels. 

The production cross section of pp —>■ dn~^ in the centre-of-mass system can be ex- 
pressed in terms of the matrix elements of the production operator as 

c?aprod _ 1 LJqEdEiE2q 1 ^ g^.j |2_ p p 

dn "(27r)2 ps ^^M^"^ >\-^,,xK,,. (4) 

Here q is the pion momentum and p the initial nucleon momentum. The energies Ei and 
E2 are the initial nucleon energies, while Uq and i^^ are the final state energies for the 
pion and deuteron, and s is the Mandelstam variable. All kinematic quantities refer to 
the centre-of-mass system. The cross section is explicitly factorized to the phase space 
factor Ppp and the sum of squared matrix elements Rpp. For the np reaction (2) one needs 
also the additional isospin factor | relative to the reaction (1) and our factorization is 
2a {np) = Pnp ■ Rnp, so that also Rnp corresponds to similar pure isospin one matrix 
elements as Rpp does. The notation is valid for both the differential and integrated cross 
sections. In the former the R is just angle dependent. 
The absorption cross section (3) 

is related to its inverse by detailed balance, i.e. by the factor ^p^ /q^- Here the integrated 
absorption cross section needs an additional factor ^ from the identity of the protons in 
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the case of positive pion absorption. The np reactions need the same factor | for the 
isospin already in the differential cross section, so the two total absorption cross sections 
should be equal for exact isospin symmetry. A similar factorization as above is possible 
also here with the same quantities R, but different phase space factors. Of course, the 
inverse of (2) is hardly experimentally feasible, but we keep it for symmetry. Often it is 
better to present the cross sections in terms of (Tabs and Fig. 1 facilitates a comparison 
in this presentation. 

Both the phase space factor (P) and the sum of the squared matrix elements (i?) 
depend on the charge via masses and different momenta. One sees easily the following 
relation for the difference of the production cross sections (both differential and total) 

Sa ^ 2a{np) - a{pp) = ^^^^^{Rnp - Rpp) + (Pnp - Ppp)^^^^^^- 

A similar equality holds for the differential absorption cross sections with just different 
phase space factors P. (For the total absorption cross section difference one should use 
Cabs(^p) ~ crabs(pp)-) Also the sum of the cross sections can be expressed exactly as 

2<7av ^ 2a{np) + a{pp) = ^"^ + ^^^ (i^„, + R,,) + (P„, - Ppp) ^'~^ , (7) 

where the latter term is only a high order correction. Neglecting it as a good approxi- 
mation one gets finally the relative change of the cross section simply as the sum of the 
relative changes of the phase space factor and of the squared matrix elements 

5a 5R 5P 

— ^-^ + -^ (8) 

^av -'''av av 

with Pav and Pav average quantities between the two reactions. 

From the above discussion it seems most reasonable to present the theoretical results 
in terms of the calculated ratio SR/Rgy with 

SR = Rnp -Rpp=J2^\<€\H''\ 0^"^ >np\' -\<i^^\H^\ 0^"^ >pp|^ (9) 

which can be used above to estimate the relative difi'erence in the cross sections. Replacing 
the average a"av by either app or (j„j, one can then get an estimate of the other cross section 
with a very small error. This prescription presumably minimizes the systematic effects of 
possible normalization deviations of the model calculation from experiments, since if the 
cross section is overestimated, likely also the difference is. This expectation will be borne 
out later in the discussion of the model dependence. Also the choice of the placement of 
the constant statistical factors either in P or P now actually becomes irrelevant in these 
relative quantities. 

Fig. 1 shows 5P/Pav and SP/P^v in different CclSGS clS db function of 77 = q/m^+. 
First, the dashed curve shows only the mass difference effect. At threshold with delicate 
cancellations of the oscillatory NN wave function, the long-ranged Galilean invariance 
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term oc {p + p') ■ {(Ji + (^2) []T0| yields a larger amplitude for the reaction (2). Although 
this is only a quarter of the dominant s-wave rescattering amplitude, this source is 
still responsible for about +2 percent units in 5i?/i?av Further, the pp momentum 
for the same pion momentum close to threshold is about two percent larger than the np 
momentum. Because of the wave function normalization as ~ jbipf + this makes the 
cross section of the reaction (1) 4 % smaller than (2). The wide and deep dip in Fig. 
1 is due to the A isobar, which increases both cross sections. However, for any given 
pion momentum the necessary pp kinetic energy is higher and, therefore, closer to the 
resonance as long as the energy is still below the A threshold. (At threshold the difference 
in the distance to the A is about 5 %.) Therefore, the reaction (1) becomes stronger 
than (2) just below the A region, while the effect in 5a is reversed above it, where the pp 
energy is further from the A. Here it is essential to remember that the distance of the 
AN mass from the two-nucleon mass is the same for both np and pp channels as discussed 
above before Eq. (4). This strong energy dependence of the cross section difference may 
be somewhat unexpected but perfectly understandable. In size the mass effect alone is, 
in fact, already comparable to today's best experimental uncertainties. 

However, in pp dn^ there is also the Coulomb force present. Although for the high- 
energy initial state this can arguably be neglected, it is important in the final state close 
to threshold. It may be noted that for larger pion energies both Coulomb effects become 
comparable, since the pp reduced mass is much larger that the pion mass. However, the 
Coulomb effect in the initial pp state is always smooth and rather small and we do not 
consider it here. The two become about equal in the neighbourhood of r/ ^ 1.2. 

In this work we consider the Coulomb effect in three different ways. First we take it 



into account in different pion partial waves by the Coulomb penetration factors [|12 

e 



CoiO ■ (10) 



,n=l 

Here the denominator (2/ + 1)!! present also in the series expansion of the spherical Bessel 
functions has been removed for a direct comparison with plane waves. As can be seen in 
e.g. the Appendix of the first of Refs. for plane wave pions some spherical Bessel 



functions ji{qr/2) enter the overlap integrals, where I is not necessarily the same as the 
pion angular momentum 1^^ relative to the deuteron. We multiply these by the above 
Ci{C,). As = arrij-edC^ / (ficq) ~ O.OOGS/r/ is rather small for any presently reported value 
of ?7, all penetration factors are practically the same as Co(0 = [27r^(e^'^^ — 1)^"'^]"'^/^ 
(mred is the ird reduced mass). The difference with the inclusion also of this effect is 
shown by the solid curve in Fig. 1. Below the A resonance the negative difference is 
partly cancelled off because of the Coulomb suppression of the dn'^ final state, while at 
threshold SR/Rg^^ approaches 2. 

Also for completeness, the dotted and dash-dot curves show the relative change of 
the phase space factors 6P/ Pav for the production and absorption reactions, respectively. 
The relative cross section difference 6a/aav can now be obtained simply by adding to 
SR/Rav the corresponding SP/Pav, as shown in Eq. H. 
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A relatively standard way of including Coulomb corrections has been to apply to the 
cross section the penetration factor Cq. Using only this gave results indistinguishable 
from the procedure described above. Also the use of Ci^ in various pion amplitudes 
directly did not deviate significantly. Another way (though still not exact) of including 
the Coulomb effect is to explicitly distort the Bessel functions to the Coulomb functions 
ji Fi. However, one should note that this has some difficulty as having different 
cases: either an emission of a charged pion with also a charged nucleon (as pn^) 
or a proton or a A"*" producing a 7r+, which interacts by the Coulomb force (possibly 
combined with s-wave rescattering) with the second nucleon. A reasonable estimate may 
be obtained by a simple replacement of the wave functions which should maximize the 
effect of this change. Only in case of serious disagreement with the earlier result is there 
reason for worry. However, the results obtained in this way were essentially the same as 
those with the simpler multiplicative penetration factors. 

A third way of looking at the Coulomb effects is to consider also the finite size of 
the charge distribution of the deuteron. For this purpose we have modified the Coulomb 
potential due to the charge distribution of the Reid deuteron wave function to 

Vi(r-,) = e^/'<iSJ*^<Ml) (11) 



Ifv - r/2| 



and used it in a Schrodinger equation to numerically solve the relevant pion wave functions 
to be used in the matrix elements. This approach neglects the nonspherical parts of the 
potential due to the D state, which would cause a coupling between different pion partial 
waves. The effect of charge extension is to very slightly increase the penetration (i.e. 
weaken the repulsion), but contrary to the expectation of Ref. the change is not 
large enough at any energy to warrant an additional calculation of the observables (at 
the level 1.5% for s waves and 0.5% for p waves). Table I shows the penetration factors 
Cq and Cf for a point-like and extended charge at a few low momenta. In the threshold 
parametrization CTprod = arj+Pr]^ these have been used in the past to scale the parameters 
a and (3. Even though the deuteron is rather extended by nuclear scale, its spatial 
variation is much faster than that of the pion wave function at very low energies. So 
basically threshold pions then feel a point charge. Further above threshold the Coulomb 
force loses importance. So either way the influence of the charge distribution is rather 
negligible to the reaction. This insensitivity to the details of the charge distribution also 
suggests that the effect of deuteron polarizability can be rather safely neglected. It may 
be noted that the present finite-size effects are only about a quarter of the results of Ref. 



0, which, integrating only to r^, missed part of the charge [0. 

Fig. 2 shows the corresponding changes in the analyzing power Ay at 90° for a 
polarized beam, which is rather representative, since its angular shape changes slowly 
and systematically. The effect of the mass difference alone is rather small, as can be 
seen from the dashed curve. The result remains indistinguishable from this, if also the 
Coulomb penetration factors are applied on the pion wave functions. This means, of 
course, that the effects are similar in different partial waves and do not change relative 
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phases. However, adding the asymptotic Coulomb phase ai^ to the amphtudes does 
change the relative phases and has a significant effect below 7] 0.8 as can be seen from 
the solid curve. To emphasize the origin of this effect the dotted curve shows still the 
difference Ay{ai^ = 0) — Ay{ai^ ^ 0) for the pure case; clearly the Coulomb phases 
oi^ dominate the change of Ay. 

It was argued above that the chosen presentation of the results as the ratio bRjRg^e 
would minimize the role of theoretical uncertainties and model dependence. We tested 
this in various ways. Using the Paris deuteron wave function instead of the Reid soft core 
one, threshold production ratio changes by about +0.6% and bAy by +0.01. In the Paris 
wave function there is a significantly softer hard core. At higher energies the changes are 
smaller. 

Next the initial state interaction was changed strongly by enhancing the NA 
transition potential by 20%. In the pion production energy region this changes the A^A^ 
phase shifts by several degrees and the cross section is changed by 10-25%. Nevertheless, 
at threshold the effects in both 5i?/i?ave and 6 Ay were totally negligible. In the A region 
the change in the ratio 5i?/i?ave was about -1% and in 6 Ay -0.005. Such a large variation 
of the transition potential is, in fact, not allowed. Namely, in Ref. a strong point was 
made that the A strength can be well fixed by the total pion production cross section 
in the A region, so there is not much uncertainty from this source. As an ultimate test 
the A was switched off and only the standard Reid soft core A^A^ potential used. At 
threshold this has an effect of reducing the cross section by a factor of | |Tl|. Even 
with this massive overall change the threshold value of the SR/Rave increased by only 
2 percent units to 9%, and then this quantity decreases monotonously from this value 
with increasing energy being about 3% at r/ = 1.8. This gives a measure of the relative 
overall importance of the A at threshold. At higher energies its role is more pronounced, 
as testified by the strong dip, which exclusively arises from approaching and passing 
the resonance. Here the A effect is about as striking as in the total cross section itself. 
However, as discussed above, there is little uncertainty in the role of the A. On the other 
hand, alterations of only the nucleon correlations have much less effect as was seen e.g. 
in Ref. [0]. 

Recently arguments have arisen about the importance of off-shell rescattering of the 
pion. Presently we consider only a monopole form factor in the whole rescattering pro- 



cess with the TiN scattering amplitude taken from the on-shell analysis of Ref. 
which assumes charge independence. While the form factor should be reasonable for the 
isovector scattering presumably mediated by mesons, it may not be sufficient for the 
isospin symmetric isoscalar channel. Varying the cut-off mass by a few hundred MeV 
had very little effect on the difference between the two reactions, while the cross section 
itself changed by tens of percent at threshold. However, at a deeper level the different 
masses of the intermediate mesons could have some charge dependent effect in the off- 
shell scattering amplitude, but such a calculation would require an actual model of vrA^ 
scattering, based perhaps on chiral perturbation theory. One would expect some effect 
also from the different pion masses in the intermediate pion propagator. However, both 
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of these effects belong to the class of dynamic isotensor forces outside the scope of this 
Letter. Work on such interactions is in progress. 

In summary, we have studied the changes in the total cross sections and analyzing 
powers Ay{90°) of the reactions pp — > d7T~^ and np — > dn^ caused by the "trivial" differ- 
ences in the masses and the Coulomb interaction. Depending on energy, these can be at 
a several percent level with varying sign, and in detailed comparisons of data from differ- 
ent reactions they should be taken into account as well as in calibrating neutron beams 
by pion production. While the Coulomb effect can be treated simply by multiplicative 
factors Co(^)^, the mass difference causes a nontrivially energy dependent effect, which 
requires an actual model calculation, the task accomphshed here. The angular distribu- 
tions change significantly less than the above observables setting an overall scale of the 
changes. Overall, our estimation of the uncertainties of the present calculation (within 
its scope to consider only the mass differences in kinematics and the Coulomb effect) 
would be about one percent unit for the 5i?/i?ave and 0.01 for 6 Ay. 
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TABLE 



Table 1: The penetration factors for a point-like charge and an extended distribution 
(labelled e). 



V 




^0 


C'f 


Cf 


0.01 


0.061 


0.061 


0.089 


0.089 


0.02 


0.286 


0.291 


0.319 


0.321 


0.04 


0.560 


0.568 


0.576 


0.579 


0.06 


0.686 


0.696 


0.695 


0.698 


0.08 


0.757 


0.768 


0.762 


0.766 


0.10 


0.802 


0.813 


0.805 


0.809 


0.15 


0.864 


0.877 


0.866 


0.871 


0.20 


0.897 


0.910 


0.898 


0.903 


0.30 


0.931 


0.943 


0.931 


0.936 



FIGURE CAPTIONS 

Fig. 1: Relative changes in the sum of the squared matrix elements (R) and phase space 
factors (P). Dashed curve: the mass difference effect in R; solid: the full result for R 
including also the Coulomb repulsion; dotted: the change in P for production; dash-dot: 
the change in P for absorption. 

Fig. 2: Differences between np dn^ and pp dir^ in the analyzing power Ay at 90°. 
Dotted: the mass difference effect alone (indistinguishable if also the penetration factors 
are included); solid: the Coulomb phase ai^ included; dashed: the resulting change in 
the pp — > dn'^ reaction, if the Coulomb phases are switched off. 
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